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Abstract 

In this paper we compare the results of the MACRO detector at 
Gran Sasso, that is the largest neutrino telescope operated before the 
year 2k, with theoretical predictions of the neutrino emission from some 
promising targets, such as blazars and GRB's. In particular we propose 
a new statistical method, justified by the maximum entropy principle, 
for assessing a model independent upper limit to the differential flux of 
neutrinos inferred from the measured up-ward going muon flux. This 
comparison confirms that the acceptance of a detector of up-ward going 
muons should be of the order of 1 km 2 , in order to challenge the present 
theoretical estimates of possible neutrino production in both galactic 
and extra-galactic objects. 



1 Introduction 

Large detectors of high energy (> 1 GeV) neutrinos, like MACRO in the 
Gran Sasso Laboratory [0 or SuperKamiokande @ have been operated in the 



1 



last decade of the past century, but no extrasolar astrophysical source has 
been detected. The occasion of this paper is related to the decommissioning 
of the MACRO detector at Gran Sasso, that is the largest neutrino telescope 
operated before the year 2k. The proposal for the detector was accepted in 
November 1984 |H|. Physics runs started in March 1989, as soon as the 
first supermodule of the detector has been built @, § The full detector 
H has been run from April, 1994 to December 31, 2000 when it has been 
decommissioned. A sample of a 1000 up-ward going muons has allowed the 
measurement of the flux of geophysical neutrinos collected up to September, 
1999 jjj, producing also interesting result ||. The purpose of this paper is 
to compare the upper limits obtained by MACRO with current models of 
neutrino emission from astrophysical objects. First in the following §|2] we 
will discuss the physics of neutrino detection in detectors of up-ward going 
muons. In particular in §2.1 we discuss the extrapolation of the differential 
neutrino cross section to the range of energies of interest for astrophysical 
neutrino detection, in §2.2 the conversion of neutrinos into muons in the 



rock and finally in § |2.3| we give the rationale for a new statistical method for 
assessing a model independent upper limit to the differential flux of neutrinos 
inferred from the measured up-ward going muon flux. In §|| we present some 
results obtained from the MACRO survey of neutrino induced up-ward going 
muons. In §3.1 we present the constraint to the neutrino emission from the 
Crab nebula, in §^2 that for a possible neutralino decay in the Galactic 
center, in §3.3 the constraint to diffuse neutrino background which can be 



inferred from the cumulative search of coincidences with blazars and in § |3.4 
we discuss the same topic for the GRB deriving an upper limit to the diffuse 
flux of neutrinos correlated with GRB . Finally in §|4| we discuss the prospect 
for detection of neutrinos in the next future by the planned detectors. 



2 Neutrino detection 

2.1 Extrapolation of neutrino cross section 

The detection process for muon neutrinos and anti neutrinos is the Charged 
Current (CC) scattering illustrated in Fig. |l| with double differential cross 
section 

d 2 ^) = G 2 F (s- m 2 N ) Af* 
dxdy 2tt (Q 2 + M 2 / ) 2 

x ^2xF^\x,Q 2 ) + (1 -y) F^\x,Q 2 ) +±L-£\ xF f\x, Q 2 ) \l) 

where y = 1 — E fl /E u is the fraction of neutrino energy lost in the Lab 
system, x = Q 2 /2mj^{E u — E„) the fraction of parton momentum carried 
by the struck parton, and Q 2 = xy(s — m 2 N ) the square of the transferred 
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Figure 1: The neutrino nucleon deep inelastic scattering in the parton model. 



four momentum. The last term is positive for z/s and negative for i/'s. 
The three functions F\,Fi and F3 are the structure functions (SF's) of the 
nucleon. In the framework of the quark parton model (QPM) the parton 
are identified with quarks, and the SF's can be expressed in terms of the 
probability density functions (PDF) of each quark types . The PDF's are 
the functions q(x,Q 2 ) where q = u, d, s, c,b,t and the relative anti-quarks, 
which assign the probability that a quark carries a momentum fraction x. 
In particular we have for an isoscalar target 

p uN(vN) = x { u + u + d + d ^ + 2s + 2c + 2h + 2t } 

x pvN{vN) _ x ^ u _- + d _ ± 2s T 2c ± 26 =F 2t} 

At very high energies the SF's F2 and 2xF%, are equal, but if the transverse 
momentum of the W is not negligible, they are related by the equation 

where R(x,Q 2 ) = o^joT is the ratio of the longitudinal to transversal cross 
section of the W. 

The structure functions F2 and xF% and the function R have been mea- 
sured experimentally up to E v ~ 500 GeV by the CCFR/NuTeV collabora- 
tion m 10, 11, 12]. To calculate the neutrino cross sections at higher ener- 



gies, as we are interested in this paper, the SF's should be extrapolated. The 
propagator term in Eq. (|) M^/{Q 2 + M^) 2 -> for Q 2 > M^, therefore 
the cross section at any energy will vanish for Q 2 > max[xy(s — m 2 N ),M' 2 v \. 
However the double differential cross section of Eq (|l|) must be integrated 
to obtain the differential cross section da/dy for a given y in the interval 
Q 2 /y(s - m 2 N ) < x < 1. If the minimum detectable energy for the muon 
in the detector is Eth the lower integration limit is x m i n w Q 2 /2m N E v . In 
Fig. [2] we have reported the area of the plane log 10 (l/ic), Q 2 corresponding 
to E v = 10 6 GeV where Eq. (||) must be integrated, compared with the area 
where the SF's have been measured. 
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Figure 2: Physical integration domain over which Eq. ([!]) must be integrated for 
a given energy of the neutrino. 
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Figure 3: Differential neutrino cross section at = 425 GeV measured by CCFR 
experiment. The fit is the differential cross section obtained in this paper applying 
CTEQ5 PDF's. 



In the same plot of Fig. we have reported the lines W 2 = Q 2 (l/x — 1) = 
4m| where mh are the masses of the heavy quarks c, b and t. Above these 
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lines the production of the heavy quarks will give an important contribution 
to the cross section. In fact it is well established that above ~ 100 GeV 
the contribution of the quark production to the total cross section is of the 
order of 20% £j]|. From the plot of Fig. § we are forced to conclude that 
the experimental knowledge of the physics underlying the extrapolation of 
the cross section is, up to now, rather incomplete. In the frame work of 
quantum cromodynamics (QCD), the gauge field theory which describes the 
strong interactions of colored quarks and gluons, the cross sections can be 
calculated perturbatively in ascending order of the strong coupling constant 
a s . In practice this corresponds to the solution of a set of partial differential 
equations dF{(x, Q 2 )/dln(Q 2 ) (the DGLAP equations) at fixed x, starting 
from the knowledge of the SF's Fi(x,Q 2 ) with i = 1,2,3 at a given lower 
value Qq (see e.g. Sterman G. et al. 1999 and references therein). This 
method has been checked, using the measured differential cross sections. 
For example we report in Fig. || the fit to the measured differential cross 



section obtained using PDF's calculated by the CTEQ5 collaboration 14], 
with evolution equations resulting from NLO analysis from a starting scale 
Qq = 1 GeV 2 . It is evident from this figure that the agreement is reasonable 
but not perfect [[l5|]. 

Nevertheless the extrapolation at higher energies of the neutrino cross 
section is still subjected to a non negligible uncertainty because, as is shown 
in Fig. ||, the domain of integration at higher energies extends to smaller 
and smaller values of x. Therefore, even if the evolution of SF's at fixed 
x to higher Q 2 can be calculated with the DGLAP equations, still we lack 
of information about the SF's at low x values. Various estimates, based on 
different physical assumption [16, 17, 18, 19, 20], as shown in Fig. |] disagree 



by factors up to 2-4 for E v > 10, 000 GeV. The experimental point has 
been obtained by the HI Collaboration pi] ]. The higher values of the cross 
section is given by extrapolation of the SF's for x < 10~ 5 based on the Regge 
theory |22|, 



In the same Fig. || we have reported the cross section of the reaction 
e~p — > vX by the HI experiment at HERA (Ahmed et al. 1994) which 
seems to be in better agreement with the lower cross section estimate. In 
the rest of this paper we will take as a baseline extrapolation of the cross 
section the more conservative one, shown as a solid line in Fig. |^, based on 
the CTEQ5 PDF's. 

2.2 Evaluation of conversion probability 

The probability P^ U (E U ) of detecting a muon with energy > E t h in the 
detector, originated by a neutrino with energy E u , is given, in the continuous 
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Figure 4: Various extrapolations of the total cross section for isoscalar target 
and equal mixture of and v^. Solid line is the baseline estimate adopted in 
the following of this paper (see text.). Finely dotted line is based on the small 
x structure functions measured at HERA, dash-dotted on the dynamical parton 
model, which is practically coincident with coarsely dotted line, dashed on the 
Regge theory for x < 10~ 5 (see text). 



slowing down approximation, by the integral 



P^v(E u ) = — I [ a + b(l-y)E v ]- 1 I ^-dydx (3) 
m N J J dxdy 



where m^v is the mass of the nucleon, o and b are the coefficient of the muon 
energy loss -dE^/dz = a(E,j) + b(E^) E^ Q. 
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Figure 5: Conversion probability for various E t h values. 



2.3 Model independent upper limits 

The differential up- ward going muon flux in the detector will be d^^/dE^ = 
P^ v d& u /dE u where d<& u /dE u is the spectral distribution of the neutrino flux. 
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However what is measured in an up-ward going muons neutrino telescope, 
like MACRO, is only the integrated muon flux 



<S>™{>E th ) = J Pr-^rdEv (4) 

E t h 

The conversion of the muon flux into a neutrino flux is subjected to the 
ambiguities arising from the lack of knowledge of the shape of the neutrino 
spectrum itself. The usual approach followed in the literature [[[], |26| is to 
give the upper limit to the neutrino fluxes assuming a spectrum d<& v jdE v oc 
E~ 2 . In this case the integrated neutrino flux can be estimated using an 
average conversion probability 

jff»(> E th ) 

<M> E th ) = i r°o F -2 P , F ( 5 ) 
E th JE th ^ a ^ 



We propose here to apply the maximum entropy principle [27, 28] to infer 
the unknown neutrino differential flux, with the less prior assumption on 
the spectral shape. In fact, even if the maximum entropy principle has 
been questioned for the lack of solid theoretical foundations, it gives very 
significant results in several fields of applications of spectral deconvolution 



[29|. According to this principle we assume that the probability distribution 
function of the energy of the neutrino that could have produced a muon 
in the MACRO detector will be uniform. In this case we have the ansatz 
Pfa, d<&„ 1E I ' dE u = const. Therefore, normalizing to the integrated neutrino 
flux of Eq. (|5|), we have derived the estimate of the upper limit to the 
neutrino flux which will be shown in the following. 



3 Results of the MACRO survey 

3.1 Neutrino emission from the Crab nebula 

Gamma ray spectrum of the Crab Nebula extends up to 25 TeV. Above 
~ 500 GeV the spectrum is a simple power law with a ~ 2.5. The total 7- 
ray luminosity is above 500 GeV L 7 ~ 1.5 x 10 26 erg cm -2 s" 1 . The detection 
of neutrino emission from the Crab could support the possibility that this 
high energy emission is due to hadronic interactions. Bednarek & Protheroe 
(1997) have analyzed the consequences of acceleration of heavy nuclei in the 
pulsar magnetosphere as a possible mechanism of energetic radiation from 
the Crab Nebula. The MACRO limit for the Crab Nebula is 

<f>C rab (> 1.2 GeV) < 2.5 x 10~ 14 cm" 2 s" 1 (90% C.L.) (6) 

In Fig. we compare the limit on the neutrino flux, with two extreme 
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Figure 6: Neutrino sky after MACRO (lower panel). In the upper panel is shown 
the complementary result of the AMANDA survey. 
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Figure 7: MACRO Upper Limit (90% C.L.) for neutrino emission from the Crab 
Nebula 



models. The present emission is determined by the trapped high energy 
proton's density in the nebula. Therefore the present emission is mainly 
due to acceleration in the early times after supernova explosion. Magnetic 
dipole losses determine the pulsar period at present, but the initial period is 
determined largely by gravitational losses and it has been probably shorter 
than 10 ms. Hence the two models consider two initial periods: 5 ms, and 
10 ms|0|. 
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3.2 Neutralino annihilation in the Galactic Center 



Intense emission from the Galactic center is predicted if cold dark matter is 
present there, as in current models of the dark galactic halo |HL |32|, [33|. In 
particular it could be expected, according to Gondolo & Silk (1999) that a 
massive black hole at the galactic center could redistribute the WIMP's into 
a cusp. The effect of this redistribution would be a strong enhancements of 
the self-annihilation rate of neutralinos. In Fig. |8]we report the predicted up- 
ward going muon flux in the two cases. After the complete run MACRO O 
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Figure 8: Predicted up-ward going muon flux from Galactic Center. The possible 
enhancement due to the redistribution of neutralinos is shown in the upper panel. 

the 90% C.L. upper limits on up-ward going muons from the Galactic Center 
is Qfj, < 0.34 x 10 -14 cm -2 s . This new limit is reported in Fig. |8[ It is clear 
from this figure that the predicted enhancement has not been observed. One 
possible explanation of this negative result has been given from a subsequent 
more careful evaluation, in a recent paper, of the dynamics of the cusp 



formation [34]. The MACRO negative observation supports the possibility 
that the accreting matter could not spiral fast enough by dynamical friction 
only. Thus within a Hubble time only a mild enhancement could have taken 
place. 

3.3 Diffuse emission from Blazars 

The second EGRET catalog of high-energy 7-ray sources |5j| contains 40 
high confidence identification of AGN and all appear to be blazars. This is 



the main reason for which several authors [36, 37] have proposed this type 



of object as potential powerful sources of high energy neutrinos. The proton 
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Figure 9: MACRO Upper Limit (90% C.L.) for the diffuse neutrino emission from 
blazars 



blazar model [38] both protons and electrons are accelerated and protons 
interact with synchrotron radiation produced by electrons photoproducing 
pions that decay into 7-rays and neutrinos. 

The up- ward going muons distribution, shown in Fig. |6| has been searched 
for angular coincidences in a 3° cone with each of the 181 blazars, listed in a 
recent catalogue with a declination 6 < 40°. No significant excess has 
been found for any individual source, with a cumulative average upper limit 
to the muon flux $ M (> 1.2GeV) < 5.44 x 1CT 16 cm^s -1 . A bias in favour 
of northern declinations is clearly present in the sample. In fact 178 out of 
the 233 objects listed in the catalog, have 5 > 0° and there are no known 
BL Lacs with 5 < —57°. We can estimate, very roughly, that the MACRO 
survey has covered the declination band —57° < 5 < 40°, corresponding to 
~ 3.2 7T sr of sky. From this figure we can estimate the U.L. to the diffuse 
up-ward going muons from blazars 



d <$>BLac 



< 0.98 x 10~ 14 cm~ 2 s -1 sr _1 



(90% C.L.) 



(7) 



Applying the maximum entropy method we obtain the upper limit to the 
diffuse neutrino flux reported in Fig. [| In the same plot we have reported 
the model independent upper bound that can be derived if the UHE cosmic 



rays are also accelerated in the same source [ 40 1 . The higher prediction is 
obtained, presuming a blazar's luminosity cosmological evolution following 
that of QSO's, which may be described as fqso = (1 + z) 3 at redshift z < 1.9 
and constant above. The lower prediction is for no evolution of the blazar's 
density. 



3.4 Search for neutrinos from GRB's 

Several mechanisms for production of intense HE and UHE neutrino and 7 
burst associated with the main pulse, in the 0.1 MeV have been proposed 



10 



in the literature [[II], [1^, [1^, It is somewhat surprising however that a 
phenomenon that was discovered in the keV-MeV region could be the herald 
of big activity in the multi-TeraVolt region. 

We have searched the catalogues BATSE 3B & 4B [45| containing 2527 
gamma ray bursts from 21 Apr. 1991 to 5 Oct. 1999 Q. They overlap 
in time with 1085 upward-going muons collected by MACRO during this 
period. In Fig. [l(] it is shown a plot of angular distance among the GRB and 
up-ward going muons in MACRO as a function of arrival time difference. 
We observe that if neutrino are massive, as suggested by observation of 
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Figure 10: Coincidences of up-ward going muons in MACRO with the 2527 gamma 
ray bursts from 21 Apr. 1991 to 5 Oct. 1999 (see text) . 

neutrino oscillations, the neutrino burst from cosmological GRB will be 
delayed respect to the observation of the 7 rays observed by the GRB itself. 
The time delay will be 

St = ^(t z -t ) (8) 

where (3 = y 1 E 2 — m 2 c 4 /E u . For E v » m u c 2 we can put (1 — j3)/j3 = 
l/2(E u /m u c 2 )~ 2 + 0(E u /m u c 2 )~ 4: . In Fig. O we show the expected delay 



of the massive neutrino signal, for GRB in the range z = 0.5 — 2 as a 
function of the Lorentz factor of the neutrino. In the same figure we have 
reported the normalized Lorentz factor distribution (for a source spectrum 
dN v /dE v oc E~ 2 ) of neutrinos that produce upward going muons in the 
MACRO detector. From this figure we see that a maximum delay of ~ 1000 s 
could be expected for very massive neutrinos (m„ = 10 eV). 

Therefore the time window to be searched for should to be expanded up 
to S3 1000 s after the detection of the GRB for the muon to be detected in 
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Figure 11: Time delay of the neutrino burst from the main pulse as a function of 
energies for different neutrino masses. 
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Figure 12: MACRO Upper Limit (90% C.L.) for the diffuse neutrino emission 
from GRB's. 



MACRO. However according to the current ideas on the production of the 
GRB explosion, it should not be possible to have the neutrino emission long 
before the GRB. Therefore we can assume that the coincidences between 
muons detected in the time window preceding the GRB should be very likely 
accidental. From the lower panel of Fig. [H] we estimate that the accidental 
coincidences among up-going muons and GRB with A6 < 20° will be 0.4 in 
the overlap period. From the upper panel of the same figure we could can- 
didate two events as possible true coincidences, one event after 39.4 s from 
4B950922 at an angular distance of 17.6° and another very horizontal event 
in coincidence with the 4B940527 inside 280 s at 14.9°. However the cumula- 
tive Poissonian probability that those two events are accidental is 7.9 x 10 -3 . 
We assume then that those events are expected accidental, and calculate that 
the 90% confidence interval for the unobserved true coincident upward-going 
muons [46] will be < 5.41. The corresponding upper limit to the muon flu- 
ence from the average GRB is $^(> 1.2 GeV) < 0.79xl0~ 9 cm~ 2 (90%C.L.) 
for neutrino masses m u < 10 eV. From this cumulative limit we can derive 
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an upper limit to the diffuse neutrino background from GRB's, taking into 
account the fact that the MACRO live time in the period from 21 Apr. 1991 
to 5 Oct. 1999 has been 4.62 y, that neutrinos could have been detected for 
GRB's with declination 5 < 40° (corresponding to 3.53 tt sr of sky) and that 
the exposure factor for BATSE in this declination range p|j is 55% we can 
convert the upper limit for the average burst to a diffuse background 

d $GRB 

-gr < 2.25 x 1(T 15 cm" 2 s" 1 sr" 1 (90% C.L.) (9) 
We obtain also in this case the upper limit to the diffuse neutrino emission 



from GRB's reported in Fig. 12. In the same figure we have reported the 
prediction of diffuse neutrino background produced by a cosmological dis- 
tribution of GRB's Also we have reported the flux expected from the 
intriguing possibility that before the jet could emerge from the dense stellar 
progenitor of the GRB, the accelerated protons could interact with the ther- 
mal X-ray photon leading to neutrinos with energy > 5 TeV @. In this 
model the intensity of diffuse neutrino background depends from the amount 
of collimation of the jet Qj et /4ir, because the ~ 1000 GRB's observed per 
year, are only a small fraction of the actual explosions. 



4 Conclusions 




Figure 13: A gallery of past, present and future HE neutrino telescopes. 

In Fig. we show a gallery of past, present and future HE neutrino 
telescopes, compared with a massive detector with 1 km 3 of volume. From 
this figure we have a direct impression of the smallness of the detectors 
which were operated in the past century, compared with the ones under 



13 



construction. It is also worth noticing that the preliminary results obtained 
by Amanda, from April to October, 1997 which are reported in the upper 
panel of Fig. [6|, include 1097 up-ward going muons. A number comparable 
to the number of events detected by MACRO in ~ 4.6 live years. This shows 
directly that the acceptance of the AMANDA-II detector is already more 
then a factor 10 greater then the acceptance of MACRO. 

In the previous §|3| we have practically shown how far was MACRO 
from effectively testing theoretical production models. From Fig l7| it can 
be crudely estimated that an acceptance ~ 10 4 MACRO ~ 1 km will be 
needed to challenge possibly the predictions for galactic neutrino sources. 
The situation could be slightly better for the statistical detection of emission 
from selected classes of objects, like for example the association of up-ward 
going muons with blazars, where a possible detection seems to be attainable 
with an acceptance ~ 0.1 km 2 . 
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